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ABSTRACT 


The  "EPSA"  (Elasto-Plastic  jShell  Analysis)  code  has  been  developed  for 
the  analysis  of  shells  in  an  acoustic  medium  subjected  to  dynamic  loadings 
which  produce  large  elasto-plastic  deformations  in  the  shell.  The  analysis 
includes  the  modeling  of  significant  internal  structures,  which  produce  hard  spots 
on  the  shell.  In  addition,  the  effects  of  ambient  pressure  are  considered.  This 

report  presents  the  theoretical  development  for  the  "EPSA"  code  and  a des- 
cription of  the  code  itself.  A users  manual  for  "EPSA"  is  planned  for  the 
future. 

The  structural  equations  of  motion  are  derived  from  the  principle  of 
virtual  work  and  descretized  over  the  shell  in  a manner  typical  of  finite 
element  procedures.  The  integration  in  time  of  the  equations  of  motion  are 
done  explicitly  via  a central  difference  scheme. 

The  nonlinear  Donnell-Vlasov  kinematic  equations  of  shell  theory  are 
used.  Plate  strain-displacement  relations  are  established  by  a two  dimensional 
finite  difference  scheme. 

Two  special  features  have  been  incorporated  into  "EPSA"  in  order  to 
obtain  a major  gain  in  the  efficiency  of  the  calculations.  First,  a self 
consistent  plasticity  theory  for  shells  has  been  developed  directly  in  terms 
of  the  stress  resultants  thereby  avoiding  conventional  " through- the-thickness" 
integrations.  Second,  a modification  of  the  basic  quadrilateral  element  has 
been  made  using  finite  difference  techniques  in  which  the  rotational  degrees 
of  freedom  are  removed  from  the  nodal  points.  As  described  in  the  report, 
both  procedures  result  in  a marked  increase  in  computational  efficiency,  parti- 
cularly for  cases  in  which  large  systems  are  to  be  analyzed. 
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The  fluid-structure  interaction  is  accounted  for  by  means  of  the 
Doubly  Asymptotic  Approximation  (DAA)  expressed  in  terms  of  orthogonal 
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I INTRODUCTION 


In  recent  years,  considerable  progress  has  been  made  in  the  development 
of  methodology  and  computer  codes  for  the  analysis  of  the  response  of  structures  in 
an  acoustic  medium  under  dynamic  loadings.  Generally,  two  types  of  problems 
have  been  studied:  (1)  the  analysis  of  structures  in  which  both  the  basic 
shell  structure  and  the  internal  components  act  elastically  under  the  dynamic 
loadings  [1],  [2];  and  (2)  the  analysis  of  structures  which  undergo  large 
elasto-plastic  deformations  under  the  dynamic  loadings. 


This  report  is  concerned  with  the  analysis  of  large  elasto-plastic 
deformations  of  stiffened  shells  under  dynamic  loadings.  Previous  reports 
by  Weidlinger  Associates  that  have  been  issued  on  this  subject  include  Ref. 
[3]  which  is  concerned  with  the  development  of  an  elasto-plastic  theory  for 
the  analysis  of  stiffened  shells  and  Ref.  [4]  which  provides  an  overview 
of  the  numerical  procedures  that  have  been  employed  on  the  dynamic  analysis 
of  elasto-plastic  shells. 

The  complexity  of  the  problems  encountered  in  developing  methodology 
for  the  elasto-plastic  analysis  of  such  shells  with  internal  structure 
requires  that  a numerical  discretization  procedure  be  utilized.  For  this 
purpose,  the  "EPSA"  (Elasto-Plastic  Shell  Analysis)  code  has  been  developed. 
"EPSA"  is  a modified  finite  element  code  which  incorporates  a number  of 
specific  features  which  are  geared  to  the  efficient  analysis  of  submerged 
shells  subjected  to  shock  Loadings. 


The  specific  objectives  which  guided  the  development  of  "EPSA”  are 
as  fallows: 

1)  Analysis  of  shells  in  an  acoustic  medium  subjected  to  both  high  and  low 


frequency  shock  loadings. 

2)  Efficient  modeling  of  elasto-plastic  behavior. 

3)  Inclusion  of  large  displacement  effects  to  analyze  dynamic 
buckling  situations  and  post-buckling  behavior. 

4)  Efficient  treatment  of  the  structure-media  interaction  problem. 

5)  Modeling  techniques  which  can  be  used  for  complex  shells  of 
arbitrary  geometry,  including  stiffeners  and  internal  structures. 

The  internal  structures  often  produce  "hard  spots"  which  materially 
affect  the  motions  of  the  shell  structure. 

In  a code  which  is  to  be  a useful  analysis/design  tool,  these  objectives 
must  be  met  in  an  efficient  fashion  to  facilitate  its  application  to  the 
solution  of  large  problems  involving  complex  structures. 

The  objective  of  this  report  is  to  present  the  basic  theory  and 
methodology  upon  which  the  EPSA  code  is  based.  The  specific  features  and 
aspects  of  "EPSA"  will  be  discussed  in  detail  in  the  Sections  which  follow. 
Two  of  these  features  have  been  utilized  to  obtain  a major  gain  in  the 
efficiency  of  the  calculations  (both  with  respect  to  accuracy  and  to  running 
time),  and  are  worthy  of  special  mention  at  this  time.  First,  an  elasto- 
plastic  shell  theory,  defined  in  terms  of  the  moments  and  direct  force 
resultants  on  the  shell,  has  been  developed.  Ref.  [3],  and  utilized  in 
EPSA.  This  is  done  in  place  of  the  usual  "through-the-thickness"  integration 
techniques  for  elasto-plastic  analysis.  The  use  of  the  elasto-plastic 
theory  results  in  a material  increase  in  the  efficiency  of  the  computational 
procedure  particularly  when  applied  to  complicated  problems. 

The  second  special  feature  involves  a modification  of  the  basic  quadri- 
lateral element  by  means  of  finite  difference  techniques,  in  which  the 
rotational  degrees  of  freedom  are  removed  from  the  nodal  points.  This 
results  in  a modified  twelve  noded,  twenty  degree  of  freedom  element  in 


which  the  degrees  of  freedom  are  all  displacements. 
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Several  advantages  are  gained  by  the  elimination  tf  the  rotational 
degrees  of  freedom.  First,  the  set  of  differential  equations  to  be  solved 
in  time  is  reduced  since  only  the  translational  degrees  of  freedom  appear 
at  each  node.  Second,  the  present  quadrilateral  element  results  in  much 
simpler  computations  than  in  any  generally  used  conforming  finite  element.  In 
addition,  one  avoids  the  possibility  of  an  ill  conditioned  mass  matrix  and 
the  related  unduly  stringent  stability  limitations  associated  with  the 
high  frequencies  produced  by  the  rotational  mass  terms. 

The  inclusion  of  the  aforementioned  special  features  makes  "EPSA" 
especially  useful  for  the  analysis  of  complicated  structures,  and  in 
particular  those  with  significant  hard  spots  introduced  by  internal  equip- 
ment. 

Comparisons  of  "EPSA"  results  with  both  analytical  and  experimental 
results  for  a series  of  problems  involving  the  dynamic  loading  of  elasto- 
plastic  shells  in  vacuo  are  presented  in  this  report. 
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1 1 ..  EQUATIONS  OF  MOTION  OF  THE  STIFFENED  SHELL 


The  equations  of  equilibrium  are  written  in  the  form  of  the  principle 
of  virtual  work: 


’ 

{s}  '{6e}  dR  - 

{p)T{<SU>  dR  + 

R 

R R 

p{U}{<5U}  dR  = 0 


(2.1) 


with. 


{U>  = (u1#  u2>  w)T 

{s}  * (nu>  n22,  n12,  m1:l,  m22,  m12) 

(e)  = (exl»  e22’  2el2 * kll*  ‘22’  2k12) 


The  scheme  used  to  solve  the  above  variational  problem  is  as  follows; 


The  surface  of  the  region  is  covered  by  a quadrilateral  me3h,  each  element  of 

which  has  an  area  A^.  The  integrals  over  R are  replaced  by  the  sums  of 

integrals  over  A^.  The  summation  over  A^  in  turn  is  performed  by  sub- 

(k) 

dividing  the  element  into  regions  A^  as  shown  in  Fig.  1.  Therefore, 

N 4 

{s}T{6e}dR=  l l }.A.(k)  = 

JR  i=l  k=l 

l l ] {(Sq}  A(k)  = {F }T{ 6q } (2.2) 
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p{U}T{6u}  dR  = 
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! I p{q)^{6<t) / (k)  - [MJ  {q}T{Sq} 
i»l  k=l 


(2.4) 


where  {q } is  the  nodal  displacement  vector  for  the  structure,  [M]  is  the 
lumped  mass  matrix,  (P)  represents  the  vector  of  external  forces  acting 
on  the  nodes  of  the  structure,  {F}  is  the  vector  of  equivalent  internal 
grid  point  forces  and  N is  the  number  of  elements. 
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The  principle  of  virtual  work  is  therefore  transformed  to  the 
following  system  of  ordinary  differential  equations: 


[M]  {q}  * I [{P)i  - {F>±] 
i*l 


(2.5) 


The  general  solution  procedure  for  Eq.  (2.5)  is  described  in  Section 
III.  An  explicit  numerical  scheme  is  utilized  for  the  forward  step 
integration  in  time. 
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III  GENERAL  SOLUTION  PROCEDURE 

The  system  of  equations  for  the  nodal  displacements  of  the 
structure,  Eq.  (2.5^  is  integrated  in  time .utilizing  the  following 

central  difference  scheme  which  is  an  explicit  method.  For  the  solution  of 
problems  involving  the  treatment  of  shock  waves  across  the  structure, 
accuracy  requirements  preclude  the  use  of  large  time  steps.  For  such 
problems  explicit  time  integration  methods  are  particularly  optimal.  The 
velocity  of  node  i at  time  step  t^+^  is  computed  by 

jt  mk>  <»•» 

where  element  forces  {Fl^  are  summed  over  all  elements  k framing  into 
node  i,  is  the  mass  of  the  node,  and  {P } ^ are  the  externally  applied 
forces  on  the  nodes. 

The  following  procedure  is  employed  at  each  time  step  of  the  solution 
phase  for  each  element  of  the  shell,* 

i)  Nodal  point  velocities  are  updated  by  the  application  of 
external  loads  occurring  within  the  time  step.  These  loads 
are  time  varying  in  either  functional  or  discrete  form. 

ii)  Boundary  conditions  are  enforced  where  applicable. 

iii)  The  strain  increment  occurring  within  the  time  step  is  computed 
from  the  nodal  point  displacements  via  the  strain-displacement 
relations.  * 

iv)  The  increment  in  the  element  stress  state  is  calculated  by  use 
of  the  constitutive  equations. 

v)  Internal  stress  state  is  converted  into  equivalent  grid  point 
forces. 

vi)  Nodal  point  velocities  are  updated  by  the  application  of  grid 
point  forces  occurring  within  the  time  step. 

Note  that  the  formulation  of  the  equations  is  in  the  initial  configuration. 


All  equations  are  solved  in  the  initial  geometry  and  all  references  are 
to  the  initial  coordinate  directions  and  initial  areas  (Lagrangiatv) . 
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IV  KINEMATIC  EQUATIONS 


The  Donnell-Vlasov  nonlinear  kinematic  equations  of  shell  theory 
are  employed.  In  orthogonal  coordinates,  the  strain-displacement  relations  read 
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The  underlined  terms  represent  geometric  nonlinearities. 
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V DISCRETIZATION 

An  arbitrarily  shaped  structure  is  divided  into  its  constitutive 
parts  called  "sheets".  Each  sheet  is  a curved  section  of  shell  with  an 
arbitrary  number  of  nodes  and  elements  (Fig.  2).  The  shape  of  the 
sheet  is  limited  to  a surface  that  can  be  described  by  a smooth  con- 
tinuous function  without  any  interior  discontinuities  in  its  slope. 

There  can  be  no  corners  or  edges  within  a sheet. 

Thus  a cylinder  with  flat  end  caps  would  consist  of  three  sheets 
a circular  cylindrical  sheet  and  a planar  sheet  for  each  end  cap  (Fig.  3). 
Three  sheets  are  required  to  specify  the  structure  because  of  the  edge 
that  occurs  between  the  cylinder  and  the  flat  end  caps. 
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Multi-sheet  capability  consists  of  assuring  the  following  boundary 
conditions  on  the  nodes  along  a shared  edge: 

1)  compatibility  of  displacements  of  nodes  along  the  edge. 

2)  compatibility  of  rotation.!  of  nodes  along  the  edge. 

3)  equilibrium  of  momenta  at  the  nodes  along  the  edge. 

Each  sheet  is  further  sub-divided  into  elements.  The  elements 
within  a sheet  can  be  of  any  arbitrary  quadrilateral  organization. 

Discrete  stiffener  elements  are  available  (Ref.  [4])  to  model  any 


stiffeners  which  exist  on  the  shell. 
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Earh  arbitrarily  shaped  quadrilateral  shell  element  is  defined  by  four 

corner  nodes,  with  each  node  having  three  translational  and  no  rotational 

*) 

degree  of  freedom  terms  . In  order  to  represent  bending  behavior  (second 
derivative  terms)  eight  nodes  not  contiguous  with  the  element  are  also 
used  (Fig.  4) . 

Each  element  accesses  twelve  nodes  and  has  twenty  degrees  of  freedom: 
three  translational  degrees  of  freedpm  for  each  of  the  four  inner  nodes, 
and  one  degree  of  freedom  (displacement  normal  to  the  surface)  for  each  of 
the  eight  exterior  nodes.  Thus,  the  nodal  displacement  vector  of  an 
element  i is 

{q>i  * (u  , u2,  u3,  u , vr  v2,  v , v , w , w , w , w , w5...w  ) (5.1) 

For  improved  accuracy,  each  element  is  further  divided  into  four 

1 

regions  for  the  computation  of  strains  in  each  region. 

The  strain  in  element  i is 


{e}±  = ({e^,  {e2>,  {e^,  {e^}^ 


where  e^  is  the  strain  in  region  1 and 


{e,}  = (e  , e , 2e  , K , K , 2K  ) 
1 xx  yy  xy  xx  yy  xy 


The  element  strains  are  related  to  the  nodal  displacements  as  follows; 


(iek}iS[Vi  {«q>i-«Bk]1+  [Bk]i  } {«q}± 


(5.2) 


(5.3) 


*)  Mnnoi  » _ll 


EPSA's"  element  formulation  differs  from  that  of  conventional  finite  element 
codes  where  rotations  and  curvatures  are  considered  as  additional  degrees 
of  freedom  and  interpolating  functions  are  used  to  establish  the  strain- 
displacement  relations. 
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[B ] T is  a function  solely  of  the  element  geometry.  [B]"  is  a function 
of  both  the  element  geometry  and  the  normal  displacements.  [Bj”  represents 
the  nonlinear  terms  in  the  strain  displacement  relations  which  are  the 
extendi*1"  terms  for  the  range  of  moderately  large  deflections.  These 
nonlinear  terms  account  for  the  finite  displacement  gradients  3w/3£^  and 
3w/3£^  the  strain-displacement  relations.  {6e}  and  {6q}  are  the 
strain  and  displacement  increments. 


The  Donnell-Vlasov  nonlinear  kinematic  equations  relating  strain  and  displacement 
increments,  written  in  Cartesian  coordinates  are: 


= 36u  + Aii  + 1”  36w 
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where  the  underlined  terms  account  for  finite  displacement  effects. 

First  derivative  (membrane)  terras  of  the  [B]  matrix  are  computed 
by  mapping  the  quadrilateral  element  via  the  Jacobian  into  a local  s-t 
coordinate  system  where  a linear  shape  function  is  assumed. 
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NT  = J [(l-s)(l-t),  'i-fsXl-t),  (l+s)(l+t),  (l-s)(l+t)]  (3.7) 


The  second  derivatives  (bending  components)  are  expressed  in  terms 


of  discrete  nodal  displacements  via  an  irregular  finite  difference  technique. 


A two  dimensional  Taylor  series  expansion  in  irregular  shaped  meshes  is 


employed.  This  technique  has  been  used  in  the  solution  of  the  large 


deflection  response  of  a flat  membrane,  Ref . [6  ] and  the  large  deflection 


response  of  a thin  shallow  spherical  shel],  Ref.  [7  ]. 


Consider  a point  P with  coordinates  x , y and  five  neighboring 
o oo 


points  P^(i  = 1>  2,... 5)  with  coordinates  x^,  y^.  The  Taylor  expansion 


of  the  normal  displacement  w re'sults  in  the  following  five  equations: 


”i  ■ "o <*i  - *o>  +!f  <>-1  - v <»i  - v2 

2 2 (5.8) 

,13  w,  N 2 , 9 w , w 

+ 2^2<yi  ' V + 3157  <xi  - *o)(yi  - V 


where  all  the  derivatives  of  w are  taken  at  P , (Fig.  4a). 


In  matrix  form. 


(a)  = [A]  {Aw} 


(5.9) 


where  Aw  = (w.  - w , - w , . . . w.  - w ) 

~ l o L o 5 o 


(5.10) 


T _ ,9w  3w  ^w  92w  92w  . 

~ 3x’  3y’  9x2’  9y2’  3x3y> 


(5.11) 


(xr  XD>*  (yi  - y0)  ••••(*!  - x0>  <yx  - y0> 


(xs  " V’  (ys  - yD>-  • • • (*5  - V(ys  - y„) 


(5.12) 


The  above  expressions  can  be  expressed  in  terms  of  the  nodal  values 


of  w as  follows: 


{a}  * [AA]  {w} 


(5.13) 


where 


{w}  * (w , w , v . . . w.) 

~ O A Z 3 





(5.14) 


-12- 


and 


[AA]  = 


“?Ali 

3 

' ] 


-1 


-1 


-IK, 


-1 


- j 


5 3 


,-1 


(5.15) 


Noting  that  the  strain  components  of  interest  consist  of  the  last 
three  components  of  the  vector  {a},  only  the  last  three  rows  of  the  matrix 
[AA]  are  needed.  With  the  pivotal  point  Pq  labeled  as  P^,  we  now  have 


{ek}i  = [AAk]  {?}i 


(5.16) 


where  w consists  of  w^  and  the  displacements  at  five  neighboring  points. 
The  selection  of  these  neighboring  points  is  not  entirely  arbitrary;  the 
following  argument  reveals  the  condition  which  must  be  imposed  on  the 
nodal  coordinates  x^,  y^  . 

The  expressions  (5.13)  and  (5.16)  yield  convergent  approximations  of 

2 

the  derivatives  3w/3x.  . . 3 w/3x3y  provided  that  the  matrix  [A]  is  non- 
singular. With  xq  = yQ  * 0,  the  singularity  of  [A]  means 


V yi  • 


Vi 


*5’  y5  * * * X5y5 


det  [A] 


(5,17) 


Thus  the -matrix  [A]  is  singular  if  there  are  four  such  real  numbers 
a,  3,  y and  6 (not  all  of  them  equal  to  zero)  such  that 

2 2 

+ ctyi  + Bxi  + yj^  +<5xiyi  = 0 (5.1 

for  i = 1,  2. . . 5. 

In  other  words,  if  the  five  neighboring  points  are  located  on  any 
conic  section  (ellipse,  hyperbola,  parabola  and  two  straight  lines)  the 
matrix  [A]  becomes  singular.  These  possibilities  rarely  occur  and  can  be 
easily  recognized  and  corrected. 

Note,  for  the  case  of  a rectangular  element,  the  derivative  expressions 
become  a staggered  finite  difference  scheme.  First  derivatives  are  com- 
puted between  nodes  and  second  derivatives  at  nodes. 

Two  options  for  choosing  elements  are  available- in  the.  EPSA  code.  An 
option  exists  to  employ  a generalized  quadrilateral  element  which  uses  the 
formulation  discussed  previously.  A second  option  exists  to  employ  a 
rectangular  element  which  uses  a staggered  finite  difference  scheme  as 
discussed  in  Ref.  [4].  A sheet  may  be  discretized  with  any  combination 

i r 

of  generalized  quadrilateral  elements  and  rectangular  elements.  * j‘ 


i 
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VI  Sh£LL  CONSTITUTIVE  EQUATIONS  - ELASTO-PLASTIC  SHELL  THEORY  DEVELOPMENT 

The  shell  constitutive  equations  relate  the  stress  resultant  rate 
vector  to  the  shell  strain  rate  vector.  This  is  formulated  in  matrix 
notation  as 

{s}  = [D]  {4}  (6.1) 

where  [D]  is  the  elasto-plastic  tangent  stiffness  matrix. 

The  shell  constitutive  equations  used  differ  from  the  classical 
elasto-plastic  theories  in  that  the  formulation  involves  shell  stress 
resultants  rather  than  stresses  at  points  throughout  the  thickness  of  the 
shell.  This  avoids  the  necessity  of  computing  and  storing  stresses  through 
the  thickness  of  the  shell  and  results  in  considerable  savings  in  computer 
storage  space  and  processing  time. 

However,  the  stress  resultants  N. . and  M, . of  shell  theory  are  not 

ij  ij 

sufficient  to  describe  the  state  of  stress.  Certain  higher-order  moments 
must  be  combined  with  the  stress  resultants  to  form  the  dynamic  variables 
of  the  problem. 

The  constitutive  relations  consist  of  a yield  condition,  a strain 
hardening  law  and  a flow  rule. 

The  stress  components  at  the  top  and  bottom  surfaces  of  a shell  are 
expressed  in  terms  of  the  stress  resultants  as 

N 6M 

a.-.  (6.2) 


ij 


h " h2 


with  the  plus  and  minus  signs  applying  to  the  top  and  bottom  fibers  of 


i 1 


the  shell. 


The  initial  yield  surface  equation  is  established  by  substituting 
the  above  relations  into  Mises  yield  condition 


72  (0112  + °222  - °11°22  + 30122)  - 1 


This  results  in 


with 


F = 1M  + XVt  + 2|X«J  " 1 

o N M NM 


XN  “ „ 2 (N11  + N22  “ N11N22  + 3N12  ) 

N 

O 

XM  = TT  (M11  + “22  M11M22  + 3M12  * 

M 

O 

INM  = N M (N11M11  + N22M22  " 2 N11N22  " 2 N22M11 
o o 


+ 3N12M12> 


where 


N = 0 h and 
o o 


0 h 

M = 

o 6 


As  the  loading  of  the  shell  is  increased  beyond  yield,  more  of  the 
wall  cross-section  plasticizes  until  eventually  limiting  values  of  the 
stress  resultants  are  reached.  This  process  is  illustrated  in  Fig.  5 
for  a simple  beam  element.  In  terms  of  the  stress  resultants,  therefore, 
the  shell  exhibits  a hardening  like  behavior  even  though  the  material  is 
modeled  as  elastic-ideally  plastic.*^ 


An  elasto-plastic  material  model  (hardening  in  the  stress-strain  curve) 
will  produce  a hardening-like  behavior  in  both  the  force  and  moment 
resultants  Fig.  6 . An  elasto-plastic  material  model  is  currently 
being  incorporated  into  "EPSA". 


fr 
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A limit  surface  is  constructed  by  assuming  a linear  combination  of 

I,  L,  and  I„,,.  Coefficients  for  these  three  terras  are  determined 
N M NM 

empirically  to  produce  a satisfactory  approximation  for  the  limit  surface. 
The  expression 

fl  ■ h * ? :h  + lm 


(6.6) 


which  represents  the  limit  condition  exactly  for  the  three  special  loading  cases 
of  (1)  membrane  forces  only,  (2)  bending  moments  only  and  (3) 

M = M22’  Ni?  = M19  = 0 was  cnosen* 


The  variable  yield  condition  which  describes  the  subsequent  yield 
surfaces  as  the  loading  path  moves  from  the  initial  yield  surface  toward 
the  limit  surface  is  generated  in  the  following  manner.  A variable  yield 
surface  of  the  form 


F = 


I«  + I * + «T  = 1 


N 


M 


NM 


(6.7) 


is  assumed  where 


M*  ■ -h  [<Mu"Mll‘>2  + (M22-M22‘)2  ' (MirMll*)  (M22-M22*> 

M 

o 


* 2 1/2 

+ 3(M12-Mi2*)  ] , 


(6.8) 


The  residual  moments  represent  "hardening  parameters"  and  are 
defined  by  the  following: 

If:  F - 1 and  "ij  > 0 

* Mo  *a 

Then:  dM,^  * 2(1  - F^)  £ 2 


ij 


o F. 


M 


If:  F<l.rj|-  + «ij  i 0 


Then  dM. . 
ij 


(6.9) 
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The  variables  F and  F„  are  defined  as 

S M 


3F 


3F 


3F 


Fs  ■ «8o^>  + <8o  §i£>  + <8o  35^> 

+ (M  + (Mo  W~^  + (Mo  3 

o 3M^  o 9^22  ° °“l2 


F„  - I<Mo  ^ + (tI0  ^ + <H0  35^>  1 


2 1/2 

] 

2 1/2 


(6.10) 


The  formulation  of  the  shell  constitutive  equations  is  completed 
with  a statement  of  the  elastic  law  and  the  flow  rule.  Noting  that  the 
shell  strain  vector  is  composed  of  an  elastic  and  plastic  portion, 

e * e'  + e" 

where  e'  is  the  elastic  strain  and  e"  is  the  plastic  strain. 


(6.11) 


The  fpllowing  elastic  law  is  assumed 

s * E(e  - e") 

where  the  elastic  matrix  E is  the  usual  she] i stiffness  matrix. 

The  plastic  strain  rates  are  defined  vVa  an  associated  flow  rule 

•„  _ , 9F 


(6.12) 


where 


C /V  rv 

- dS 

# ss 


(6.13) 


,3F.  _.3F.  . 3F  . . .3F. 

(3^  E(3^>  ’ (3i*>  A(3l> 


and 


M Fc 

A - 2(1  - F ) ^ J*_ 

° fm 


A more  detailed  discussion  of  the  constitutive  theory,  including 
numerical  results, is  presented  in  Ref.  [3j. 


' i 
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VII  STRUCTURE-MEDIUM  INTERFACE 

With  the  view  of  establishing  in  "EPSA"  the  capability  of  modeling 
complicated  structures  containing  internal  components,  an  approxiamtion 
of  the  structure-fluid  surface  interaction  which  uncouples  the  acoustic  medium 
from  the  structure  was  desirable.  Such  an  uncoupling  scheme  would  allow 
most  of  the  computer  capability  to  be  used  in  a realistic  modeling  of  the 
shell  and  its  internal  configuration. 

The  uncoupling  procedure  used  is  the  Doubly  Asymptotic  Approximation  (DAA) 
developed  by  Geers  (Ref.  [8])  and  by  Mnev  and  Pertsev  (Ref.  [9]),  Ranlet 
et  al,  Ref.  [2])  have  employed  the  DAA  method  with  a natural  modal  expansion 
technique  to  predict  accurately  the  elastic  response  of  ring-stiffened  cylinders 
with  internal  equipment  subjected  to  shock  loading.  The  DAA  therefore  has 
been  used  to  analyze  elastic  small  deflection  phenomena. 

The  inclusion  of  the  DAA  into  "EPSA"  extends  the  theory  into  the  large 
defleccion,  elasto-plastic  range. 

The  nature  of  many  transient  structure-media  interaction  problems  is 
a rapidly  applied  load  followed  by  a low-frequency  response.  The  DAA 
produces  exact  results  in  the  low  and  high  frequency  ranges.  For  short 
times,  DAA  reduces  to  a plane  pressure  wave  approximation  and  for  long 
times  it  reduces  to  a virtual  mass  approximation. 

The  DAA  imparts  upon  the  structural  model  surface  loading  composed  of 
incident  and  radiated  waves.  The  form  of  'this  loading  is 
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M Q. 


?i  ■ Vi + Ai  i rw 

k=l  Tt 


(7.1) 


where 

til 

i)  F is  force  imparted  on  ic  node  by  the  fluid. 

ii)  is  incident  pressure  obtained  by  an  empirical  representation 
of  the  explosive  loading. 

iii)  is  a coefficient  equal  to  the  total  wet  area  of  the  shell. 

til 

iv)  ijt  (s  ) is  the  k Surface  Expansion  Function  (S.E.F.) 

^ *_t. 

evaluated  at  the  i node  (total  of  M S.E.F.'s). 

v)  is  the  generalized  force  in  the  normal  direction  for  the  kfc^  S.E.F. 

The  generalized  fluid  forces  are  expanded  along  the  surface  of  the 
shell  by  means  of  Surface  Expansion  Functions  (S.E.F.).  This  technique 
has  been  successfully  employed  by  Ranlet  et  al  (Ref.  [2]).  A direct 
integration  technique  for  the  fluid  force  is  an  alternative  means  of 
solution  which  is  being  considered  for  the  extension  of  the  theory  to 
arbitrary  geometries,  Ref.  [10].  The  Surface  Expansion  Functions  represent 
an  orthogonal  set  of  functions  over  the  shell  used  to  circumvent  the 
poorly  conditioned  behavior  of  the  DAA  equations  when  using  expansions  in. 


terms  of  normal  mode  components  (Ref. [11]). 


The  generalized  fluid  force  is  obtained  from  a system  of  coupled 
first-oder  differential  equations  written  as; 

\ lvxk  - i.  ‘V1 1 V‘ " t \ i(v  c)  WdA } <7,: 

J53**-  n & a 


I r 
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where 


vik  ■ 1 1 Vv  WdA 
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V^.  is  the  incident  velocity  transmitted  by  the  fluid  upon  the  shell. 


c is  the  sound  velocity  in  fluid.  Y^j  is  a virtual  mass  coefficient  obtained 


from  the  solution  of  an  incompressible  steady-state  problem  in  which  the 
normal  displacement  of  each  surface  expansion  function  is 
applied  on  a cavity  having  the  same  shape  as  the  structure. 


Equation  (7.2)  is  solved  in  time  with  a first-order  integration  scheme 
as  follows: 


{Q}  1 = (pfcyk)  [{V}  1 - [y]  {Y}  1 - [¥]  {q}  *] 


(7.4) 


with 


{Y}  1+1  . {Y}  1 + At  {Q}  1 


(7.5) 


where  the  following  arrays  are  defined  as: 


i)  Generalized  Fluid  Force  Vector 


{Q}  = (Yl 


ii)  Incident  Velocity  Vector 

(v} ' t j,  W WAi 


i-1 


iii)  Surface  Expansion  Function  Array 


m = 


y,, 


The  generalized  fluid  force  is  applied  to  the  structural  nodes  as  a 


; 1 


! \ 


i \ 


dynamic  loading. 
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vm  COMPUTATIONAL  FEATURES 

The  need  to  realistically  model  complex  structures  containing  internal 
equipment  has  stipulated  that  efficiency  in  computation  capability  be  a 
premium.  Therefore,  the  development  of  EPSA  has  proceeded  with  the  goal 
of  being  an  efficient  analytical  tool.  To  this  end,  basic  premises 
and  theories  have  been  introduced.  These  include: 

1)  An  explicit  integration  in  time  scheme,  eliminating  the 

need  for  assembly  and  inversion  of  large  matrices  at  each  time  step. 

2)  A "Simple"  quadrilateral  element  is  used.  This  introduces  a 
large  number  of  mass  points  in  the  structure  to  achieve  an 
accurate  solution  for  shock  and  wave  propagation  problems. 

3)  The  elimination  of  rotational  degrees  of  freedom  which  cause  stringent 
stability  limitations  due  to  the  high  frequencies  produced  by  the 
rotational  mass  terms. 

4)  An  elasto-plastic  theory  employing  stress  resultants,  thereby 
eliminating  the  need  for  through  the  thickness  integration. 

With  the  incorporation  of  these  features  in  EPSA  an  efficient  means 
of  analysis  has  been  achieved. 

Coding  and  computations  have  been  carried  out  on  two  versions  of  EPSA. 

A small  core  version  which  is  run  on  the  CDC  6600  machine  and  a large  core 
version  which  is  used  on  the  CDC  7600  machine. 

For  the  small  core  version  (CDC  6600)  computational  time  is  five 
milliseconds  per  time  step  per  element,  A "typical"  problem  involving  1000 
elements  and  500  time  steps  will  consist  of  a run  time  of  fourty  minutes. 

For  each  element  of  the  structures  thirty  words  of  data  are  stored,  thereby 
limited  the  maximum  number  of  elements  per  structure  to  2500. 

For  the  large  core  version  (CDC  7600)  computational  time  is  one 
millisecond  per  time  step  per  element.  The  above  "typical"  problem  there- 
fore consists  of  a run  time  of  ten  minutes.  There  is  effectively  no  limit 
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on  the  number  of  elements  per  structure. 

A re-start  provision  has  been  incorporated  into  EPSA  thereby  providing 
the  user  with  an  effective  means  of  checking  intermediate  results. 


IX  IMPLEMENTATION 


During  all  phases  of  EPSA  development,  code  verification  was  enacted 
to  provide  a level  of  confidence  and  supportability  in  the  code.  Subsequent 
to  the  inclusion  of  each  new  feature  in  the  code’s  capability,  calculations 
were  made  to  test  the  code  against  analytical,  experimental  and  numerical 
results. 

Reference  [3]  contained  checks  of  the  elasto-plastic  portion  of  the 
code.  Included  were  comparisons  of  moment-curvature  relations  using  the 
EPSA  shell  theory  versus  the  classical  through-the-thickness  technique. 
Reference  [4]  reported  on  a circular  cylindrical  shell  subjected  to  small 
elasto-plastic  deformation.  EPSA  Analysis  compared  favorably  with  similar 
computations  using  the  DYNAPLAS  code,  Ref.  [16], 

Analyses  were  made  by  EPSA  to  verify  the  large  deflection  capability 
for  both  elastic  and  elasto-plastic  stress  states  due  to  dynamic  loading. 

Some  such  results  will  be  presented  as  follows: 

The  large-deflection  static  response  of  a fixed-ended  beam  subjected 
to  a point  load  was  investigated.  Ref. [12],  Results  were  obtained 
assuming  (1)  linear  behavior  and  (2)  nonlinear  behavior  (large  deflections). 
Problem  specifics  and  results  are  shown  in  Fig.  7. 

Excellent  agreement  was  obtained  with  analytical  results.  The 
computer  model  consisted  of  6 elements  per  half  length  of  beam.  Computational 
time  was  less  than  two  minutes  per  calculation. 

The  instability  characteristics  of  a fixed-ended  arch  was  also 
investigated,  Ref.  [13).  The  magnitude  of  a concentrated  load  was  varied 
to  determine  what  the  ultimate  load  capacity  of  the  structure  is.  Problem 
specifics  and  results  are  shown  in  Fig.  8.  A comparison  with  a Finite 
Element  Solution  by  Marcal,  Ref.  [14],  is  also  shown.  Both  results  are  in 
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agreement.  The  computer  model  consisted  of  6 elements  per  half  length  of 
arch.  Computational  time  was  less  than  two  minutes  per  calculation. 

The  most  worthy  determination  of  whether  "F.PSA"  can  indeed  provide 
accurate  and  reliable  predictions  for  large  deflection,  elasto-pj ash ic 
structural  response  is  determined  by  comparisons  of  EPSA  predictions  against 
relevant  experimental  data.  To  this  end,  the  elasto-plastic  lar^  deformation 
transient  and  permanent  response  of  a circumferentially  stiffened  cylindrical 
panel  was  investigated. 

R.W.H.  Wu  and  E.A.  Witmer  at  M.I.T.  performed  experiments  on  a 
integrally  stiffened  clamped-edge  6061-T6  aluminum  cylinder  panel  subjected 
to  impulsive  loading  by  the  sheet  explosive  loadiftg  technique.  Ref.  [15]. 

The  geometric  properties  of  the  panel  are  shown  in  Fig.  9.  To  insure 
ideally  clamped  edge  conditions,  the  panel  was  machined  out  of  an  aluminum 
block  leaving  a rectangular  collar  for  fixity.  For  repeatability  and 
confirmation  of  results  three  such  panels  were  tested.  A sheet  of  high 
explosives  centered  on  the  upper  surface  of  the  panel  was  detonated 
providing  the  impulsive  loading. 

Analysis  of  the  structure  by  "EPSA"  consisted  of  a 225  element  mesh 
configuration  for  the  quarter  model.  The  nominal  yield  of  46  ksi  was 
increased  to  60  ksi  to  account  for  strain  rate  effects.  Initial  conditions 
consist  of  initial  nodal  radial  velocites.  These  are  obtained  by  equating 
the  impulse  imparted  by  the  detonated  sheet  to  the  total  impulse  experienced 
by  the  structure.  The  time-step  size  used  in  the  integration  phase  was 
0.5  psec.  Total  computer,  run  time  was  less  than  15  minutes. 
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Figure  10  shows  comparisons  for  the  deformed  panel  shape  of  EPSA 
versus  experimental  values.  Comparisons  of  strain  time  histories  are 
seen  in  Fig.  11.  Excellent  agreement  was  obtained  for  both  deformations 


and  strains. 
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X CONCLUSIONS 

A practical  method  for  predicting  the  inelastic  response  of  structures 
in  an  acoustic  medium  under  dynamic  loadings  is  presented  in  this  report. 

The  "EPSA"  code  has  been  developed  to  serve  as  an  efficient  tool  which  is 
particularly  aimed  at  meeting  the  complexities  involved  in  modeling  such 
structures  with  major  internal  components.  EPSA  development  has  proceeded  with 
computational  efficiency  as  the  major  objective. 

The  excellent  agreement  whcih  has  been  achieved  by  EPSA  in  comparison 
with  experimental,  analytical  and  computational  results  has  ver4Cied  the 
correctness  of  the  approaches  employed  in  the  code, 

With  this  developing  level  of  confidence  in  the  EPSA  code,  additional 
capabilities  are  currently  being  incorporated.  At  various  stages  of 
development,  these  are: 

(1)  An  improved  multi-sheet  capability  to  provide  for  the  analysis 

of  cylindrical  shells  with  hemispherical  or  conical  end  closures. 

(2)  Inclusion  of  ambient  pressure  acting  on  the  3hell. 

(3)  Implementation  of  a general  initial  stress  and  deformation  state 
in  the  structural  model, 

(4)  Provision  for  out-of-plane  bending  of  stiffeners  (tripping). 

(5)  A consideration,  in  the  analysis,  of  major  internal  structures, 
etc. 

In  conjunction  with  the  expansion  of  EPSA's  capability,  additional 
experimental  work  is  being  planned.  Yarious  scaled  and  configured  models 
will  be  subjected  to  dynamic  loadings  leading  to  large  elasto-plastic 
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deformations.  The  behavior  of  these  increasingly  complex  models  will 
provide  a rational  basis  for  the  verification  of  EPSA  as  additional 
capabilities  are  provided  in  the  code. 
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